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ABSTRACT 

A quantitative spectral analysis of 24 A supergiants in the Sculptor Group spiral 
galaxy NGC 300 at a distance of 1.9 Mpc is presented. A new method is introduced 
to analyze low resolution (~ 5A) spectra, which yields metallicities accurate to 0.2 
dex including the uncertainties arising from the errors in T e g (5%) and log g (0.2 dex). 
For the first time the stellar metallicity gradient based on elements such as titanium 



Based on VLT observations for ESO Large Programme 171.D-0004. 
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and iron in a galaxy beyond the Local Group is investigated. Solar metallicity is mea- 
sured in the center and 0.3 solar in the outskirts and a logarithmic gradient of —0.08 
dex/kpc. An average reddening of E(B-V) ~ 0.12 mag is obtained, however with a 
large variation from 0.07 to 0.24 mag. We also determine stellar radii, luminosities 
and masses and discuss the evolutionary status. Finally, the observed relationship be- 
tween absolute bolometric magnitudes M& / and flux weighted gravities gp = g/T e t is 
investigated. At high temperatures the strengths of the Balmer lines depends solely on 
the flux-weighted gravity, which allows a precise direct determination of log gf with 
an accuracy of 0.05 to 0.1 dex. We find a tight relationship between M^ / and log gF 
in agreement with stellar evolution theory. Combining these new results with previous 
work on Local Group galaxies we obtain a new flux weighted gravity luminosity rela- 
tionship (FGLR), which is very well defined and appears to be an excellent alternative 
tool to determine distances to galaxies. 

Subject headings: galaxies: — distances galaxies: abundances — galaxies: stellar 
content — galaxies: individual (NGC300) — stars: early-type 



Introduction 



It has long been the dream of stellar astronomers to study individual stellar objects in distant 
galaxies to obtain detailed spectroscopic information about the star formation history and chemo- 
dynamical evolution of galaxies and to determine accurate distances based on the determination 
of stellar parameters and interstellar reddening and extinction. The most promising objects for 
such studies are massive stars in a mass range between 15 to 40 M & in the short-lived evolutionary 
phase, when they leave the hydrogen main-sequence and cross the HRD in a few thousand years 
as blue supergiants of late B and early A spectral type. Because of the strongly reduced absolute 
value of bolometric correction when evolving towards smaller temperature these objects increase 
their brightness in visual light and become the optically brightest stars in the universe with abso- 
lute visual magnitudes up to My = —9.5 rivaling with the integrated light brightness of globular 
clusters and dwarf spheroidal galaxies. 

There has been a long history of quantitative sp ectrosc o pic st udies of these extreme objects. 
In a pioneering and comprehensive paper on Deneb iGrothl (|196lh was the fi rst to obtain s tellar 
param e ters and deta iled chemical composition. This work was continued by IWolfl (|197l|) . IWolf 
(119721). I W olf (1973) in studies of A supergiants in the Milky Way and the Magellanic Clouds. 
Kudritzkil (| 19731) using newly developed NLTE model atmospheres found that at the low grav- 
ities and the correspondingly low electron densities of these objects effects of de partu r es from 
LTE can become extremely important. With strongly improved model atmospheres IVennl (|1995a|) . 
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Vennl (I1995bl) and lA ufdenbe rg et al.l (I2002|) continued these studies in the Milky Way. Most re- 
cently ,|ftzybillaeLaL ( 2006 ) used ve ry detailed NLTE line formation calculations including ten 



thousands of lines in NLTE (see also [Pr zvbilla et al.l l2000l IPrzybilla et al.ll2001al IPrzybilla et al 



2001bl IPrzybilla & Butlerll200 U IPrzybilla Il2002|) to determine stellar parameters and abundances 
with hitherto unkown precision (7^ to < 2%, log g to ~ 0.05 dex, individual metal abundances to 
~ 0.05 dex). At the same time, utelizing the power of the new 8m to 10m class tele scopes, high 
resolutio n studies of A s u perg iants in ma ny Local Group galax ies were carried out 



Venn et al. 



(SMC), lM cCarthy_etaLl (119951) (M33), lM cCarthy_etaL| d 19971) ( M31), 



by 



Venn et al 



feoOll) (NGC 6822), Ivenn et all J2003h (WLM). and iKaufer et"aD (|2004j) (Sextans A) 



Venn1(ll999l) 
d2000b (M31), 



yielding invaluable information about the stellar chemical composition in these galaxies. 

The concept to go beyond the Local Group and to study A supergiants by means of quantitative 
spectr osco py in galax i es out to the Virgo cluster has been fi rst presented bylKudritzki. Lennon & Puis 



(1995) and Kudritzki ( 1998 ). Following-up on this idea, Bresolin et al. ( 200l|) and iBresolin et al. 
(|2002l) used the VLT and FORS at 5 A resolution for a first investigation of blue supergiants in 
NGC 3621 (6.7 Mpc) and NGC 300 (1.9 Mpc). They were able to demonstrate that for these dis- 
tances and at this resolution spectra of sufficient S/N c an be obtained allowing for the quantit ative 
determination of stellar parameters and metallicities. iKudritzki. Bresolin & Przvbilla d2003l) ex- 
tended this work and showed that stellar gravities and temperatures determined from the spectral 
analysis can be used to determine distances to galaxies by using the correlation between absolute 
bolometric magnitude and flux weighted gravity gf = g/T^ (FGLR). However, while these were 
encouraging steps towards the use of A supergiants as quantitative diagnostic tools of galaxies 
beyond the Local Group, the work presented in these papers had still a fundamental deficiency. 
At the low resolution of 5 A it is not possible to use ionization equilibria for the determination of 
Tgff in the same way as in the high resolution work mentioned in the previous paragraph. Instead, 
spectral types were determined and a simple spectral type - temperature relationship as obtained 
for the Milky Way was used to determine effective temperatures and then gravities and metallici- 
ties. Since the spectral type - T e « relationship must depend on metallicity (and also gravities), the 



method becomes inconsistent as soon as the metallicity is significantly different from solar (or the 
gravities a re larger than for lumino sity class la) and may lead to inaccurate stellar parameters. As 
shown by Evans & Howarthl (|2003|) . the uncertainties introduced in this way could be significant 
and would make it impossible to use the FGLR for distance determinations. In addition, the metal- 
licities derived might be unreliable. This posed a serious problem for the the low resolution study 
of A supergiants in distant galaxies. 

In this paper we present a solution to this problem and provide the first self-consistent de- 
termination of stellar parameters and metallicities for A supergiants in galaxies beyond the Local 
Group based on the quantitative model atmosphere analysis of low resolution spectra. We apply 
our new method on 24 supergiants of spectral type B8 to A5 in the Scultor Group spiral galaxy 
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NGC 300 and obtain temperatures, gravities, metallicities, radii, luminosities and masses. Note 
that while we have included a number of late type B supergiants in this work, we refer to the whole 
sample as "A supergiants " in the following, for r easons of simplicity. We will compare some of our 



results with the work by lUrbaneja et al.l (|2005|) . who performed a quantitative study of six early 
B-type supergiants (BO. 5 to B3). We will refer to those objects as "B supergiants". The reason for 
this simplifying destinction is the difference in the analysis method. The temperatures of the B su- 
pergiants (in this nomenclature) are determined from the ionization equilibria of helium and silicon 
and metallicities result from the fit of individual spectral lines of a-elements such as oxygen, mag- 
nesium, silicon and carbon, whereas the temperatures of the A supergiants are obtained from the 
Balmer jump and the metallicities are determined from the spectral synthesis of many overlapping 
lines of different species including heavy elements such as iron, chromium and titanium. 

The paper is organized as follows. After a brief description of the observations in section 
2 we provide a detailed discussion of the analysis method, i.e. the determination of T e ff, log g, 
metallicity, E(B-V), radius, luminosity and mass from the low resolution spectra and HST pho- 
tometry, in section 3. Section 4 presents the results of the analysis together with a discussion of 
their uncertainties. The metallicity of NGC 300 obtained from this work is discussed in section 
5 and compared with previous work. Section 6 investigates stellar properties and evolutionary 
time scales. The new FGLR is presented in section 7 together with a detailed discussion of the 
uncertainties. Conclusions and outlook towards future work are given in section 8. 



2. Observations 



The spectroscopic observations were obtained with FORS1 at the ESO VLT in multiobject 
spectroscopy mode in two consecutive nights of excellent seeing conditions (typical 0.7 arcsec, 
but long spell s of 0.4 to 0.5 arcsec) on 2000 September 25 and 26 as part of the Auracaria Project 
(|Gieren et al.l (|2005b() ). Four fields (A, B, C, D) were observed in five exposures, each lasting 
45 minutes. The spectral resolution is approximately 5 A and the spectral coverage extends from 
shortward of the Balm er jump to H/3 or bey ond for most of our targets. The observational data set 
has been described by iBresolin et al.l (I2002|) . The paper also contains finding charts, coordinates, 
photometry (ESO/MPI 2.2m WFI) and spectral types. Of the 62 confirmed supergiants found 
with spectral types ranging from late O to F, 29 were of spectral type B8 to A5 - the spectral 
type of interest for this study - and were not severely contaminated by HII emission or composite 
(indicating a contribution of a star of different spectral type). 24 of those had a S/N per pixel 
of about 40 or better, sufficient for our quantitative analysis. This is the sample listed in Table 
1 and selected for th e work presented in the following. The target designation is the same as in 
Bresolin et al.l (|2002[) . Note that a quantitative spectral analysis of the early type B supergiants 



(spectral type BO to B3) has already been published by lUrbaneja et al 



(2005). 



The analysis method presented in the next section makes use of the strength of the observed 
Balmer jump D B . (Fortunately, FORS was developed with high UV throughput as a design goal, 
which turns out to be an enormous advantage for our study). To measure Dg we selected one 45 
minute exposure for each of the observed fields (usually the exposure with the best seeing) and 
used standard stars to obtain a flux calibration of our spectra, where special care was given to the 
calibration in the near UV. While the single exposures are somewhat noisy, they do allow for a 
rather accurate measurement of Dg (see next section). Slit losses due to the 1 arcsec slitlets can be 
corrected for by renormalizing the observed flux calibrated spectrum to the HST B-band photome- 
try. They are unimportant to first order, since Dg is a relative measurement of fluxes shortward and 
longward of 3647 A. However, we have to worry about second order effects caused by the wave- 
length dependence of seeing and resulting wavelength dependent slit losses. Fortunately, FORS is 
equipped with an atmospheric dispersion compensator, which was used for our observations, thus, 
already correcting for slit losses caused by atmospheric dispersion. The effects of wavelength de- 
pendent seeing were addressed in the following way. We measured the FWHM of our spectra on 
the detector perpendicular to the direction of dispersion as a function of wavelength. This provided 
us with an estimate of the seeing as a function of wavelength. Assuming a Moffat-function (with 
exponent /3 = 2) for the PSF we were then able to calculate slit-losses at every wavelength from 
the measured seeing FWHM using the simple formulae 



old n + (a/w) 2 ) - 5 
Q.1169FWHM(X) 



pnew _ p 



a 



(1) 



w is the slit half-width corresponding to 0.5 arcsec. 

The relative flux corrections obtained in this way are small over the wavelength range, which 
we use for measuring or fitting the Balmer jump Dg (3585 A to 3876 A). The logarithm of ratio 
of the correction factors at these two wavelength are 0.016, 0.070, 0.024 and 0.004 for fields A, 
B, C, D respectively, smaller or comparable to the stochastic errors caused by the photon noise. 
Assuming a Moffat-function with /3 = 3 yields very similar results. 

The original photometry for our targets was obtai ned with the ESQ 2. 2 m and the Wid e Field 
Imager as a result of a multiepoch survey for Cepheids (|Bresolin et al ]l2002h . bresolin et all J2004h 
have used these data for a study of blue supergiant photometric variability. No significant variabil- 
ity has been found for our targets. 

Almost all of our targets were observed with HST/ ACS in an inve stigation of the effects 
of blending on the Cepheid distance to NGC 300 (Br esolin et al.1 120051) . which were found to 
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be negligible. With A supergiants being significantly brighter than Cepheids it is obvious that 
blending and multiplicity should not be an issue for this stu dy, as is also confirmed by a comparison 
of the ground-based and HST/ACS V and I photometry. iBresolin et al.l (|2005l) detected an off- 
set of 0.1 mag between ESO/WFI and HST/ACS B-band photometry, which they attributed to a 
calibration uncertainty of the ground-based data. We will use the HST/ACS photometry for our 
targets, where available. For three objects, we do not have HST data. In these cases, we use the 
ESO/WFI B, V, I data but apply a 0.1 mag correction to the B-band. 



3. Analysis Method 

For the quantitative analysis of the spectra we will use the same com bination of line blank eted 
model atmospheres and very detailed NLTE line formation calculations as Przybilla et al. ( 2006|) in 



their high signal-to-noise and high spectral resolution study of galactic A- supergiants. The model 
atmosphere assumptions and the individual atomic models used for the NLTE calculations have 
been described in this paper, which also contains a variety of tests with regard to the accuracy and 
consistency that can be obtained. In the spectral range to be used in our work (3500 A to 5500 A) 
these calculations reproduce the observed normalized spectra and the spect ral energy distribution , 



including the Balmer jump, extremely well. For the four objects studied, iPrzybilla et al.l (|2006) 
calculated a relatively small grid of models centered around the expected parameters for each star 
followed by more line formation calculations to obtain element abundances, once the stellar param- 
eters were determined. The large sample of 24 objects in our study requires a different approach. 
We calculate an extensive, comprehensive and dense grid of model atmospheres and NLTE line 
formation covering the potential full parameter range of all the objects in gravity, effective temper- 
ature and metallicity. We then use the synthetic spectra and energy distributions obtained in this 
multi-dimensional parameter space to constrain the stellar parameters. The gravities and tempera- 
tures of the grid are displayed in Fig. [Q 

At each grid point we have a set of models with line formation calculations for the following 
metallicities [Z] = log Z/Z : -1.30, -1.00, -0.85, -0.70, -0.60, -0.50, -0.40, -0.30, -0.15, 0.00, 0.15, 
0.30. Z/Z is the metallicity relative to the sun in a sense that the abundance for each element 
is scaled by the same fa ctor relative to its solar abundance. Solar abundances were taken from 



Grevesse & Sauvall (|1998|) . We note that by scaling abundances in this way we are not able to 
detect deviations from the solar abundance pattern but only general changes in metallicity. We are 
confident that it will be possible to extent this work to at least a study of the relative ratio of a- to 
iron-group elements as a function of galactocentric distance. But at this stage we restrict our work 
to a discussion of stellar metallicities and the metallicity gradient and postpone a more detailed 
study to follow-up work. 
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A determination of metallicities always requires an estimate of the photospheric mictrotur- 
bulence velocity y t normally obtained from the constraint that weak and strong lines of each ion 
should yield the same element abundance. At our spectral resolution it is difficult to constrain 
mictroturbulence in this way, thus, we rely on the results of the high resolution work done for A 
supergiants in the Milky Way and Local Group gala xies, which, for t unately, shows a very con 



sistent pattern of v f as a function of stellar gravit ies (Przybilla 



1995a, 



Venn 1995bL IVenn et alJl200d I Venn et all 12001 



200 1 IPrzvbilla et al 



Venn et al. 



2003, Kaufer etal 



2006, Venn 



2004). At 



low gravities the values of v f are about 8 km/sec and decrease to lower values of about 4 km/sec at 
higher gravities. The values of \ t adopted for this study, accordingly, are also displayed in Fig.Q] 

The analysis of the each of the 24 targets proceeds in three steps. First, the stellar param- 
eters (Tgff and log g) are determined together with interstellar reddening and extinction, then the 
metallicity is determined and finally, assuming a distance to NGC 300, we obtain stellar radii, lu- 
minosities and masses. For the first step, a well established method to obtain the stellar parameters 
of supergiants of late B to early A spectral type is to use ionization equilibria of weak metal lines 
(OI/II; Mgl/II; NI/II etc.) for the determination of effective temperature T e g and the Balmer lines 
for the gravities log g. The most recent high resolution work cited in section Q] has very convinc- 
ingly shown that very accurate results can be obtained in this way, when detailed non-LTE line 
formation calculations are used. However, at the low resolution of 5 A the weak spectral lines of 
the neutral species disappear in the noise of the spectra and an alternative technique is required to 
obtain temperature information. 

One simple way is to use the information from stronger lines in the spectrum which define the 
spectral type and then to apply a spectral type - temperature ca libration based on the high resolution 



quanti tativ e spectroscopic work. This techn i que w as used by [Bresolin et al.l (|200l|) . iBresolin et al. 



(120021) and iKudritzki. Bresolin & Przybillal (I2003|) . The disadvanteage of this method is obvious. 
As soon as the metallicity of the objects investigated with this low resolution method is signif- 
icantly different from the metallicity of the objects which are used to calibrate the spectral type 
effective temperature relationship, substantial errors can be introduced. For instance, objects at a 
lower metallicity will have weaker characteristic metal lines and will thus be assigned a spectral 
type earlier than an object of the same effective temperature but higher metallicity. Thus, with- 
out a priori knowledge o f the metallicity this method might indeed introduce large uncertainties. 
Evans & Howarthl (12003b have investigated this effect for A supergiant of luminosity class II in the 
Magellanic Clouds by studying the Call lines and found significant effects. However, our targets 
are more luminous (luminosity class I) and the spectral type classification used in the papers men- 
tioned above did not use the Call lines, but the metal line spectrum in the range from 4400 to 4700 
A. Thus, an extended investigation is needed. 



Fortunately, the comprehensive model atmosphere and line formation grid described above 
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makes it easy to investigate the degeneracy of the metal line spectrum with temperature and metal- 
licity. All we have to do is to follow the iscontours of constant Balmer line equivalent widths in 
the (log g, log Tgff) plane (see Fig. [5]) (along which the Balmer lines as a gravity indicators have 
identical profiles at 5 A resolution) and to compare synthetic metal line spectra of models with 
different T e f and [Z]. A rather dramatic example obtained in this way is shown in Fig. [2] with metal 
line spectra of three models at different pairs of T e g and log g, where the Balmer line profiles are 
exactly the same. The models have different metallicity [Z]. The hottest model at T e g = 10500 K 
with a metallicity slightly higher than solar ([Z] = 0. 15) has an almost identical metal line spectrum 
as the coolest model at T e g = 8750 K with a metallicity significantly lower than solar ([Z] = -0.50). 

Thus, while all three models have exactly the same Balmer line profiles and would fit the 
observed Balmer lines of a hypothetical star to be analyzed, we would have a hard time to inde- 
pendently determine temperature and metallicity from the rest of the spectrum, which looks very 
similar for all three models. Fortunately, the strength of Hel A 4471 allows, in principle, to dis- 
tinguish between the hottest model and the two cooler ones, and this will be one of the alternative 
techniques that we will use in the following. However, it seems impossible to extract the temper- 
ature difference between the two cooler models from the spectra displayed. We note the slightly 
enhanced strengths of the Till lines at low temperature, but using these lines as an independent 
temperature indicator while simultaneously fitting the Balmer lines would be based on an assump- 
tion on the abundance ratio of titanium to iron, and, thus, the ratio of a -elements to iron-group 
elements. While the investigation of the possible variation of the ratio as a function of galaxy evo- 
lution, is beyond the scope of this paper (see above), it will certainly become important in future 
work. We, thus, need a temperature indicator independent from any a priori assumptions about 
general metallicity and relative abundance ratios of metals. This is in particular important, as the 
spiral galaxies, which we intend to investigate with our spectroscopic technique, are assumed to 
have significant metallicity gradients and possibly also gradients, for instance, in their a to iron- 
group element ratios. 

A very obvious way out of this dilemma is the use of the spectral energy distributions (SEDs) 
and here, in particular of the Balmer jump Dg- Fig. [3] and @] show, how the degeneracy of the 
temperature vs. metallicity diagnostic is broken. Very obviously, the SEDs and Dg's of the three 
models with almost identical metal line spectra are significantly different. While the observed 
photometry from B-band to I-band will be used to constrain the interstellar reddening, we can use 
Dg as a temperature diagnostic. For stars of spectral type A this idea is certainly not new and goes 
back to the classical work by Chalonge and Divan i n the 1950's. F or A supergian t s it has been 



successfully investigated and applied by iGrothl (| 196 II) . I Ay din I (| 1972|) . and lKudritzkil (| 19731) . 



The basis for our diagnostic technique is given in Fig. [5] It shows isocontours of constant H5 
equivalent width and of constant values of the Balmer jump D B . The Balmer jump is defined as 
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below 

D B = (log(F l ™ 8 ))-(log(Fl h ° rt )) 

(log(F l ° n8 )) = UogF ° + logF ° + logF ° + logF °}/4 

\ 6\ x J/ 6 3?82A 6 3gl4A 5 3847A 5 3876A J/ (2) 

(log(F[ hort )) = {^log(F Xi )}/N, 3585A < h < 3627 'A 

i=\ 

The wavelengths for the individual fluxes on the longward side are selected to avoid the wings 
of the Balmer lines and stronger metal line features. In this definition, the dependence of Dg on 
[Z] is very weak and can be neglected, as an inspection of our model grid shows. The same is 
true for the strengths of the Balmer lines and, thus, the diagnostic diagram of Fig. [5]is metallicity 
independent. While the two sets of isocontours are not perfectly orthogonal, they certainly lead to 
well defined intersection points except very close to the Eddington-limit. 



3.1. Effective temperature, Gravity and Reddening 

Guided by Fig.|5]we apply the following diagnostic technique to determine effective temper- 
ature and gravity. We use the normalized spectra of our targets to interactively fit the Balmer lines 
profiles Hio to H y with our model spectra grid in the (log g,log T e g) plane (the fit curves obtained 
in this way follow very closely the isocontours of Balmer lines as in the example of Fig. [5]). Then, 
we fit the observed SED and the Balmer jump in the same diagnostic plane. Our multi-color pho- 
tometry from B-band to I-band compared with the model fluxes allows to determine reddening and 
extinction and the calibrated FORS fluxes around the Balmer jump (see section |2]) fitted with the 
model fluxes yield a second fit curve for Dg (again very similar to the isocontours in Fig. [5]). The 
intersection of the two fit curves defines temperature and gravity. 

An example is given in Figs. [6l [7J and [8] for object No. 21 of our sample. At a fixed 
value of effective temperature the Balmer lines allow for a determination of log g within 0.05 dex 
uncertainty. The uncertainty of Balmer jump Dg measured from the SED fit can be (conservatively) 
estimated to be 0.02 dex corresponding to a temperature uncertainty at the fixed gravity of 

1.55 at Tgff = 10000 K. (The values of Dg obtained from SED fit are given in Table 1). Applying 
the same fit procedure at different effective temperatures in the (log g, log T e g) plane leads to the 
fit diagram of Fig. [8j While the uncertainties at a fixed pair of T e g and log g are rather small, the 
fact that the two sets of fit curves for the Balmer lines and the Balmer jump are not orthogonal 
leads to the significantly larger error box, which is obtained from the intersection of the dashed 
curves. The uncertainty of our low resolution analysis obtained in this way are for T e ff and 

In 17 dex for log g, respectively. This is roughly a factor of two larger than what can be obtained 
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with high quality, high r esolution spectra and the use of ionization equilibria from weak lines 
(ses IPrzybilla et al.l I2006Q but still good enough to determine metallicities [Z] with a reasonable 
accuracy, as we will demonstrate below. The fit of T e g and log g also yields the value of interstellar 
reddening E(B-V) = 0.13 m ag from a compariso n of computed model fluxes with fluxes obtained 
from HST photometry (see iBresolin et al.ll2004l for the photomet r ic dat a). As for the calculation 



of the reddening law we use the formulae given by lCardelli et al.l (| 19891) and Ry = 3.1 for the the 
ratio of visual extinction Ay and reddening E(B-V). 

The accurate determination of T e g and log g is crucial for the use of A supergiants as distance 
indicators using the relationship between absolute bolometric magnitude M\, \ and flux weighted 
gravity log gp defined as 



log gF = logg- AlogT effA (3) 



where T e g A = Ziyi0000K(see lKudritzki. Bresolin & Przybillal 120031) . The relatively large un 



certainties obtained with our fit method may casts doubts whether log gf can be obtained accu- 
rately enough. Fortunately, the non-orthogonal behaviour of the fit curves in Fig.[8]leads to errors 
in T e ff and log g, which are correlated in a way that reduces the uncertainties of log gf. This is 
demonstrated in Fig.|9l which shows the corresponding fit curves of the Balmer lines and Dg in the 
(log g.log T e ff) plane. We obtain much smaller uncertainties of log gp for target No. 21, namely 
— 08 • ^ e discuss the physical reason behind this in section [6] 

We have applied this diagnostic technique on all targets, for which we were able to obtain 
observations of the Balmer jump. In two cases (target No. and 10), the set-up of the FORS 
slitlets did not allow to obtain spectral information shortward of the Balmer edge. For those, we 
use the strong temperature dependence of the observed Hel lines to constrain T e g (see Fig. [2]). For 
targets No. 2 and 4 we have used the observed Hel lines as an additional temperature criterium, 
since for both of them the intersections of the Dg and Balmer line fit curves are not well defined. 
The errors in T e g obtained in these four cases are comparable to those obtained with the Balmer 
jump for most of the other stars. The Hel lines provide a very useful alternative diagnostic tool, 
however this technique works only at high enough temperatures and low gravities. The results for 
T e ff, log g, log gf and E(B-V) for all targets are summarized in Table 1 . 

Generally, the uncertainties for T e g, log g, and log gF are similar to those obtained for our 
example target. However, for the effective temperature errors AT e ff/T e ff there is a slight trend with 
temperature as shown by Fig. [T0l which comes from the fact that the Balmer jump becomes less 
temperature dependent at higher temperatures (see Fig. [5]). On average, temperatures are accurate 
to 5 percent. This is an uncertainty rough ly twice as large as f or the analysis of high resolution and 



high S/N spectra of similar objects (see IPrzybilla et al.ll2006|) . The uncertainties for gravities and 
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flux weighted gravities are 0.2 dex and 0. 1 dex on average, again twice as large as for the analysis 
of high quality spectra. This is a result of the larger uncertainty of the temperature determination. 



3.2. Metallicity 

Knowing the stellar atmospheric parameters T e g and log g we are able to determine stellar 
metallicities by fitting the metal lines with our comprehensive grid of line formation calculations. 
The fit procedure proceeds in several steps. First, we define spectral windows, for which a good 
definition of the continuum is possible and which are relatively undisturbed by flaws in the spec- 
trum (for instance caused by cosmic events) or interstellar emission and absorption. A typical 
spectral window used for all of our targets is the wavelength interval 4497 A < X < 4607 A. 
Fig. [TTI shows the synthetic spectrum calculated for the atmopsheric parameters of target No. 21 
(the example of the previous subsection) and for all the metallicities of our grid ranging from -1 .30 
< [Z] < 0.30. It is very obvious that the strengths of the metal line features are a strong function 
of metallicity. In Fig. [[2] we show the observed spectrum of target No. 21 in this spectral region 
overplotted by the synthetic spectrum for each metallicity. We use separate plots for each metallic- 
ity, because the optimal relative normalization of the observed and calculated spectra is obviously 
metallicity dependent. We address this problem by renormalizing the observed spectrum for each 
metallicity so that the synthetic spectrum always intersects the observations at the same value at 
the two edges of the spectral window (indicated by the dashed vertical lines). The scaling of the 
observed spectrum between the two edges is done by a simple linear interpolation. The next step 
is a pixel-by-pixel comparison of calculated and normalized observed fluxes for each metallicity. 

n pix 

(0-C) 2 ([Z]) ; = [£(F([Z])f -F{[Z])f c f]/n pix (4) 
The index i refers to spectral window selected. If we define 

X ([Z})i = S/N*(0-C)([Z})i, (5) 

where S/N is the signal-to-noise ration per resolution element, we can look for the minimum 
#([Z]) ; - as a function of [Z] to determine metallicity. This is done in Fig. \T3\ which yields [Z] = 
-0.40. Application of the same method on different spectral windows provides additional indepen- 
dent information on metallicity as shown in Fig. [14] and \\5\ The average metallicity obtained from 
all windows is [Z] = -0.39 with a very small dispersion of only 0.02 dex. 

While it seems that our method yields a very accurate metallicity for a fixed pair of stellar 
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parameters T e g and log g, we also need to consider the effect of the stellar parameter uncertainties 
on the metallicity determination. This is done by applying the same correlation method for [Z] for 
models at the extremes of the error box for T e ff and log g, as given by Fig. [8] and Table 1 . This 
increases the uncertainty of [Z] to ± 0.15 dex. 

We apply the same method to all the objects of our sample. The results are given in Table 
1. In some case the uncertainties are larger than for target No. 21, which is caused by larger 
uncertainties of the atmospheric parameters or more noise in the spectra or both. 



3.3. Radii, Luminosities and Masses 

The fit of the observed photometric fluxes with the model atmosphere fluxes was used to de- 
termine interstellar reddening E(B-V) and extinction Ay = 3.1 E(B-V). Simultaneously, the fit also 
yields the stellar angular diameter, which provides the stellar radius, if a distance is adopted. Al- 
ternatively, for the stellar parameters (T e ff, log g, [Z]) determined through the spectral analysis the 
model atmospheres also yield bolometric corrections BC, which we use to determine bolometric 
magnitudes. These bolometric magnitudes then also allow us to compute radii. The radii deter- 
mined with these two different methods agree withi n a few percent. We have chosen the latter 



method for the radii given in Table 2. iGieren et al.1 (12005 a|) in their multi-wavelength study of 



a large sample of Cepheids in NGC 300 including the near-IR have determined a new distance 
modulus m-M = 26.37 mag, which corresponds to a distance of 1.88 Mpc. We have adopted these 
values to obtain the radii and absolute magnitudes in Table 2. 

In the temperature range of 8300 K < 12000 K of our model grid (see Fig. Q]) the individual 
BC-values of the models are well represented by an analytical fit formula 

BC(T eff , log g, [Z]) = f(log g F ) - 4.3 log T 4 (l + log T 4 ) +D Z (6) 
where T 4 = 7^/10000 K and 

D Z = 0.09[Z](1+0.26[Z]) 
f(log g F ) = /mflx(l - a exp{-{log gF-x min )/h)) 

and a = l-fmJfmax, fmin = "0.39, fmax = -0.265, x min = 1.075, h=0.17. The fit is accurate to 
0.02 mag over the whole range of model parameters. It allows us to estimate the uncertainties of 
bolometric magnitudes arising from the uncertainties of the stellar parameters. The errors for for 
M ho i given in Table 2 include these uncertainties together with the photometric errors. 
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We note that in most cases fit formulae provided for bolometric corrections do not include 
the effects of gravity. In the case of A supergiants this would lead to significant error residuals of 
the fit of the order of 0.4 mag. Thus, the inclusion of the function f(log gp) is essential for the 
expression of BC. 

In our model atmosphere approach to calculate bolometric corrections we convert the model 
flux at the V-band effective wavelength into a V-magnitude in the usual way using a constant, which 
transforms the flux of Vega into the Vega V-magnitude. We then use the integral over the total 
model spectral energy distribution and the effective temperature of the Sun (5777K) to calculate 
the bolometric magnitude (assuming M? , = 4.74 mag). The difference between the bolometric 
and the V-magnitudes calculated in this way then yields the bolometric correction. This approach 
provides the correct bolometric correction for Vega, which is exactly in the effective temperature 
range of the objects we are dealing with. We have, therefore, decided not to re-normalize our 
BC values to avo i d sma ll positive values at low effective temperatures, as has been suggested by 
Buser & Kurucz I (| 19781) (but see also the full discussion on this matter in this paper). As a result, 
our BC-values are ab out 0.12 mag larger than the ones given by IPrzybilla et al.l (120061) . who have 
used the approach by [Buser & Kurucz I (|1978|) . 



There are two ways to determine stellar masses. We can use the stellar gravities together 
with the radii to directly calculate masses. We call the masses derived in this way spectroscopic 
masses. Alternatively, masses can be estimated by compar ing the location of our target stars in 
the HRD with evolutionary tracks (|Meynet & Maederl 120051) . Masses obtained in such a way are 
called evolutionary masses. Both values are given in Table 2 for each target star. The uncertainties 
are discussed in section 4.2. 



4. Results 



In this section we discuss the main results of our quantitative spectral analysis. 



4.1. Interstellar Reddening 

Our quantitative spectroscopic method yields interstellar reddening and extinction as a by- 
product of the analysis process. For objects embedded in the dusty disk of a star forming spiral 
galaxy we expect a wide range of interstellar reddening E(B-V). Indeed, we find a range from E(B- 
V) = 0.07 mag up to 0.24 mag. Fig.[l6]shows the distribution of interstellar reddening among the 24 
stars of our sample. The individual reddening values are significantly l arger than the value of 0.03 



mag adopted in the HST distance scale key project study of Cepheids by LFreedman et al.l(|2001 ) and 
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demonstrates the need for a reliable reddening determination for stellar distance indicators, at least 
as long the study is restricted to optical wavelengths. The average o yer our samp l e is (E (B — V)) 
= 0.12 mag in close agreement with the value of 0.1 mag found by iGieren et al.l (|2005a|) in their 
optical to near-IR study of Cepheids in NGC 300. While Cepheids have somewhat lower masses 
than the A supergiants of our study and are consequently somewhat older, they nonetheless belong 
to the same population and are found at similar sites. Thus, we expect them to be affected by 
interstellar reddening in the same way as A supergiants. 



4.2. Stellar properties 



The original motivation for the development of our analysis technique using the Balmer jump 
was the degeneracy between the effective temperature derived from the spectral type and metallic - 
ity. In Fig. [FT] and [T8l we compare the temperatures obtained with the two different methods. We 
use the compilation by lKudritzki. Bresolin & Przybillal (|2003l) for the relationship between spectral 
type and effective temperature. While there is a clear correlation between the two temperatures, 
there is definitely a shift to higher temperatures for ^{spectral type), which is explained by the fact 
that all of our targets have a metallicity smaller than in the solar neighborhood. Fig. [T9l indicates 
that the effect becomes larger with decreasing metallicity, exactly as expected (for the discussion 
of metallicities we refer to the next section). 



Fig. [20] shows th e location of our targ ets in the (log g, log T e g) plane compared to t he early B- 

super giants studied by lUrbaneja et al.l(|2005|) . The comparison with evolutionary tracks (|Meynet & Maeder 
2005b gives a first indication of the stellar masses in a range from 10 M Q to 40 M . Three targets 
have obviously higher masses than the rest of the sam ple and seem to be on a similar evolutionary 
track as the objects studied by lUrbaneja et al.l ((2005). 



An alternative way to discuss the evolutionary status of our targets is the classical HRD. 
This is done in Fig. [2T] The evolutionary information obtained from this diagram compared with 
Fig. [20] appears to be consistent. The B-supergiants seem to be more massive than most of the A 
supergiants. The same three A supergiants apparently more massive than the rest in Fig. [20] are 
also the most luminous objects in Fig.[2T] This confirms that quantitative spectroscopy is -at least 
qualitatively - capable to retrieve th e information about absolute luminosities. Note that the fact 
that all the B supergiants studied by lUrbaneja et al.1 (120051) a re more massive is si mply a selection 
effect of the V magnitude limited spectroscopic survey by iBresolin et al.l (|2002|) . At similar V 
magnitude as the A supergiants those objects have higher bolometric corrections because of their 
higher effective temperatures and are, therefore, more luminous and massive. 



The luminosity uncertainties for our targets range from 0.02 dex to 0.06 dex (the latter rather 
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large errors appear only in a few rare cases, see Table 2). Correspondingly, Fig. [2T] allows to 
determine stellar masses based on evolutio nary tracks with a much b etter resolution than Fig. [20j 
We have used the evolutionary tracks by iMeynet & Maeden (12005) (and references therein) to 
construct mass luminosity relationships in the A supergiant stage of the form 



logL/L Q = (a- bx)x + c, x = logM/M Q - logM min /M Q (8) 

The coefficients for different evolutionary models (with and without rotation, Milky Way and 
SMC metallicity) are given in Table 4. The fit reproduces model masses for a given luminosity (in 
the range —6.5 > M\, \ > —9.7) to better than 5% or 0.02 dex. With the bolometric magnitudes 
in Table 2 we can then determine evolutionary masses, which are also given in Table 2. We have 
chosen the mass-luminosity relationship for tracks with rotation and SMC abundance, the other 
relationships give similar results. The reason why we selected Milky Way and SMC metallicities 
as a representative grid of evolutionary tracks is that only for those a complete set of tracks is 
available. 

With the small errors in luminosity, we obtain masses which formally have very small uncer- 
tainties of at most 5%. However, the masses are, of course, affected by the systematic uncertainties 
of the evolutionary tracks through the effects of mass-loss, rotation, mixing and opacity. One par- 
ticular problem may arise through blue loops from the red giant branch (RGB) or crossing to the 
blue from the RGB. Fortunately, for the mass and luminosity range of interest (12M < M < 30M Q , 
—6.5 < M bo j < —9.3) we find that the published tracks (in particular at lower metallicity) do not 
show significant loops or second crossings. 

An alternative way to determine masses is, of course, to use directly the stellar gravities and 
radii determined through the fitting of stellar spectra and SEDs. These spectroscopic masses are 
also given in Table 2. Note that we have not included the effects of centrifugal forces on the 
effective stellar grav ities to derive spectroscopic masses. While this is important for O-stars (see 



Repolust et al.1120 04). the rotational velocities of A supergiants (~ 30 km/s) are too small and the 
radii too large to cause an effect. Fig.[22]and [23]compare the masses derived by the two different 
methods. Given the fact that uncertainties of the spectroscopic masses are roughly as large as the 
gravity uncertainties (0.1 to 0.2 dex) the comparison is generally encouraging. At the high mass 
end, we have agreement within the uncertainties. At the low mass end, however, the evolutionary 
masses seem to be somewhat higher on average. There is the indication of a trend of this effect 
to become stronger with decreasing spectroscopic mass. At this point, we have no explanation for 
this result. If it is not an artefact of the spectroscopic technique applied, then it may indicate that 
the effects of blue loops are more important than the evolutionary models available to us predict. 
Blue loops increase the luminosity at a given mass ( at the low mass end the amount of mass-loss is 
still small at the RGB), thus, their neglect would lead to an overestimate of the evolutionary mass. 
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The most extreme object in the left part of Figures [22] and [23] ist target No. 17 of Tables 1 
and 2. It is the object with the lowest luminosity in the HRD (Fig. I2TT) corresponding to a mass of 
1 1 M . Its position in the (log g, log T e jf) plane of Fig. [20]is right on the track with an initial mass 
of 15 M , but with the low gravity determined through spectroscopy its radius is too small to yield 
mass high enough in comparison with the evolutionary mass. 



4.3. Metallicity Gradient 

The metallicities of the individual targets together with their galacto-centric distance are given 
in Table 1 and allow us to discuss the stellar metallicity and the metallicity gradient in NGC 300. 
Fig. [24] shows the stellar metallicities [Z] as function of galactocentric distance, expressed in terms 
of the isophotal radius, p/po- Despite the scatter caused by the metallicity uncertainties of the 
individual stars the metallicity gradient of the young disk population in NGC 300 is very clearly 
visible. A linear regression (excluding target No. 1, the outlier at p/po = 0.81 and [Z] = -0.04) for 
the combined A- and B-supergiant sample yields 

[Z] = -0.06 ±0.09 -0.44 ±0.12 p/po (9) 
for the angular gradient. On a linear scale (d in kpc) we obtain 



[Z] = -0.06 ± 0.09 - 0.083 ± 0.022 d. 



(10) 



We obtain a very similar regression, if we restrict the sample to A-supergiants only, which 
indicates that there is no systematic abundan ce difference betwee n these two groups of stars. Our 
new result confirms the result obtained by lUrbaneja et al.l (|2005f) . but with a much larger sample 
the metallicity gradient and the zero point are now much better defined. Note that the metallicities 
of th e B supergiants refer to oxygen only with a value of log[N(0)/N(H)] = -3.31 adopted for the 
sun ( Allende Prieto. Lambert & Asplundll2001|) . On the other hand, the A supergiant metallicities 
reflect the abundances of a variety of heavy elements such as Ti, Fe, Cr, Si, S, and Mg. 



With this result we can resume the discussion started by Urbaneia et al. (2005) and compare 



with oxygen abunda nces obtained from Hll-r egion emission lines. lUrbaneja et al 



(2005) used line 



fluxes published by iDeharveng et al.l (11988) and applied various different published strong line 
method calibrations to determine nebular oxygen abundances, which could then be used to obtain 
the similar regressions as above. The results for zero points and gradients c ompared to our work a re 
given in Table 3 (for a discussion of the individual calibrations we refer to lUrbaneja et al]|2005f) . 
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The different strong line method calibrations lead to significant d ifferences in the centra l 
met allicity as wel l as in the abundance gradient. The calibrations by iDopita & Evans] (I1986Q 
and IZaritsky et al.1 (119941) predict a metallicity significantly supersolar in the center of NGC 300 
contrary to the other calibrations. Our work yields a central metallicity sli ghtly smaller than solar i n 
good agreemen t with De nicolo et al.l (120021) and marginally agreeing with lKobulnicky et al.1 (|1999l) . 
Pilyuginl (|200ll) . and lPettini & Pagell (|2004|) . At the isophotal radius, 5 .3 kpc away from the center 
of NGC 300, we obtain an average metallicity significant ly smaller than solar \Z ] = -0.50, close 



to the a verage metallicity i n the SMC. The calibrations by IDopita & Evans] (119861) . IZaritsky et al. 



(I1994|) . Kobulnicky et al.l (|1999l) do not reach these small values for oxygen in the HII regions 
either because their central metallicity values are too high or the metallicity are gradients too 
shallow. 

In light of the substantial range of metallicies obtained from HII region emission lines using 
different strong line method calibrations it seems to be extremely valuable to have an independent 
method using young stars. It will be very important to compare our results with advanced work 
on HII regions, which will avoid strong line methods, and will use direct infomation about nebular 
electron temperatures and densities (Bresolin et al., 2008, in preparation). 



As mentioned before, there is a suspicious outlier in Fig. [241 We have checked the stellar 
parameters and metallicity of this object several times and have not found a reason to discard the 
high metallicity value found, which is 0.38 dex above the average value at this galacto-centric 
distance. A similar example is target No. 14 at the largest galacto-centric distance. While we 
cannot exclude that something went wrong in the analysis, we also want to state very firmly that 
the expectation of homogeneous azimuthal metallicity in patchy star forming spiral galaxies seems 
naive to us. The continuation of this type of work on other galaxies will show whether cases like 
this are c ommon or not. W e note that the high resolution spectral analysis of A supergiants in 
WLM by lVenn et al.l (|2003l) has detected a similar outlier case. 



So far, we have only tried to determine an average metallicity [Z] for each star. Together with 
the challenge of an accurate determination of stellar parameters we felt that this was a good first 
reasonable step to obtain quantitative information from our low resolution spectra. While [Z] in 
the way how our method works is mostly determined by an overlap of Fe, Ti and Cr lines, we 
have also observed that certain spectral windows are dominated by a majority of lines of specific 
element species. In future work we will try to use this observation to determine relative fractions of 
abundances of a- elements such as Mg and Ti to iron-group elements such as Fe and Cr. It will be 
an important goal to investigate whether this ratio is constant or not as a function of galactocentric 
distance. 
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5. Stellar Properties and Evolutionary Lifetimes 

The progenitor stars of the A supergiants of our s ample were obviously O-stars in a mass range 



between 15 to 25 M . The HST/ACS photometry of Bresolin et al.l (120051) allows to estimate how 



many of such progenitor stars are present in exactly the same field. Fig.[25]shows the corresponding 
combined color-magnitude diagram. Adopting our mean reddening of E(B-V) = 0.12 together 
with the corresponding mean visual extinction of Ay = 0.37, a distance modulus of 26.37 mag, 
an absolute magnitude of My = - 2.7 mag for a 15 M star on the zero age main sequence, we 
can assume that all stars brighter than my = 24.0 mag are possible progenitors of A supergiants. 
Fig.[26]shows the number of star counts in our fields in bins of 0.1 mag as a function of magnitude. 
The curve is well fitted by 



log n = a+my + 



(ID 



and a-k = 0.5622, b+ = -9.921 (see also iBresolin et 19961 bresolin et all 1998b . While the 
number of star counts starts to become incomplete at my = 21.5 mag we can extrapolate the fit 
to calculate the total cumulative number N of stars down to my = 24.0. We obtain N= 2.9 x 10 4 
for the present number of progenitor stars. Assuming continuous star formation and having 24 
confirmed A supergiants in the fields the ratio 24/2.9 x 10 4 = 8.3 x 10~ 4 provides a lower limit 
for the evolutionary lifetime in the A supergiant stage in units of the main sequence lifetime. 

The total number of A supergi ants in the obs e rved fields in the mass range between 15 to 
25 M Q can then be estimated from Bresolin et alj (|2002l) . who found a total of 167 supergiant 
candidates in the combined color magnitude diagram. They obtained spectra for 70 of those and 
detected 29 A supergiants, i.e. a fraction of 41%. Assuming the same fraction for the original 
total of 167 candidates leaves us with an estimate of 69 A supergiants. The observed ratio of 
evolutionary lifetime in the A supergiant stage to the main sequence lifetime is then 69/2.9 x 10 4 = 
2.4 x 10~ 3 . 



The corresponding lifetime ratio for evolutionary tr acks is a strong functio n of metallicity. 
Comparing calculations including the effects of rotation (IMeynet & Maeden 120051) we obtain 2 x 
10~ 2 for SMC metallicity and 2 — 8 x 10~ 4 for Milky Way metallicity. With metallicities for most 
of our target stars between these extremes we conclude that our empirical lifetime estimate agrees 
with the prediction of the evolutionary models. 



- 19- 



6. Flux Weighted Gravity - Luminosity Relationship (FGLR) 

As discussed in the previous sections, massive stars with masses in the range from 12 M to 
40 M evolve through the B and A supergiant stage at roughly constant luminosity (see Fig. |2TI) . 
In addition, since the evolutionary timescale is very short when crossing through the B and A 
supergiant domain, the amount of mass lost in this stage is small. This means that the evolution 
proceeds at constant mass and constant luminosity. This has a very simple, but very important 
consequence for the relationship of gravity and effective temperature along each evolutionary track. 
From 

L oc R 2 T* S = const. ;M = const. (12) 

follows immediately that 

Moc g R z ocL(g/T* s ) =Lg F = const. (13) 

Thus, along the evolution through the B and A supergiant domain the "flux-weighted gravity" 
gf = glT^ should remain constant. This means each evolutionary track of different luminosity in 
this domain is characterized by a specific value of gf. This value is determined by the relationship 
between stellar mass and luminosity, which in a first approximation is a power law 

LocM x . (14) 

We note from section 4.2 that the mass-luminosity relationship has higher order terms, but for 
simplicity we use this as a first order approximation to obtain 

L l ~ x ~ (g/T? s ) x . (15) 
With the definition of bolometric magnitude M^ \ oc — 2. 5 log L we then derive 



-M bo i = a FGL R{log gf - 1.5)) + bfGLR ■ (16) 

This is the "flux-weighted gravity - luminosity relationship" (FGLR) of blue supergiants. 
Note that the proportionality constant afGLR is given by the exponent of the mass - luminosity 
power law through 



aFGLR = 2.5x/(l -x) , 



(17) 
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for instance, for x = 3, we obtain cipglr = —3.75. We have chosen the zero point of the 
relationship at a flux weighted gravity of 1.5, which is in the middle of the range encountered for 
blue supergiant stars. 

We can use the mass-luminosity relationships of different evolutionary tracks (with and with- 
out rotation, for Milky Way and SMC metallicity) as discussed in section 4.2 and characterized 
by the coefficients given in Table 4 to calculate the different stellar evolution FGLRs, which are 
displayed in Fig. [27] Very interestingly, while different evolutionary model types yield somewhat 
different FGLRs, the differences are rather small. The difference is largest between the two sets 
of stellar evolution models for Milky Way metallicity with and without the effects of rotation. 
The models with SMC metallicity and rotation yield a FGLR rather close to the case with Milky 
way metallicity and rotation. Only at higher luminosity, where the dependence of mass-loss on 
metallicity starts to become important, is a larger difference from the Milky Way FGLR predicted. 
This is also the parameter domain, where the curvature of the FGLR becomes more significant, as 
to be expected from the non-linear mass-luminosity relationship of section 4.2. The set of stellar 
evolution models at SMC metallicity but without the effects of rotation yields a FGLR very close 
to the one with rotation and the same metallicity. 



Kud ritzki. Bresolin & Przybillal (|2003|) were the first to realize that the FGLR has a very in- 
teresting potential as a purely spectroscopic distance indicator, as it relates two spectroscopically 
well defined quantitities, effective temperature and gravity, to the absolute magnitude. Compiling a 
large data set of spectroscopic high resolution studies of A supergiants in the Local Group and with 
an approximate analysis of low resolution data of a few targets in galaxies beyond the Local Group 
(see discussion in previous chapters) they were able to prove the existence of an observational 
FGLR rather similar to the theoretically predicted one. 

With our improved analysis technique of low resolution spectra of A supergiants and with 
the much larger sample studied in our new approach we can now resume the investigation of the 
FGLR. We will do this in three steps. We will first discuss the determination of flux weighted 
gravities log gF from the Balmer lines and will then introduce the new F GLR obtained for NGC 



300. In this second step we will also include the B supergiants analyzed by lUrbaneja et al.1 (12005). 
Finally, we will compare the new results obtained for NGC 300 with the previous work on Local 
Group A supergiants to find out whether a universal FGLR exists. 



6.1. Balmer lines and Flux weighted Gravity 



In Fig. [28] we show the dependence of the strength of the Balmer line H§ characterized by its 
equivalent width W^(H5) on gravity log g and effective temperature T e g. Obviously, W^(H5) de- 
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pends on both, gravity and temperature, and the non-horizontal isocontours displayed in Fig. [5] are 
the result of this dependence. The situation changes dramatically, if one looks at the dependence 
°f W^(H5) on flux- weighted gravity log gF as shown in Fig. [29] Now all the curves for higher 
effective temperatures are on top of each other and form a temperature independent relationship, 
which leads to practically horizontal isocontours in the {log g, log T e g) plane Fig. [30] All the other 
Balmer lines from Hy to Hio show a very similar behaviour. 

This has important consequences for the diagnostics of flux-weighted gravities. For objects 
hot enough the uncertainties of the effective temperature determination do not affect the determi- 
nation of effective gravities. A simple measurement of the strengths of Balmer lines either through 
line profile fitting or by simply measuring the equivalent widths gives log gp with rather high 
precision, as we have already pointed out in section 3. 

In terms of distance determination using the FGLR this is very encouraging. It means that er- 
rors in T e ff stemming from the application of an inaccurate relationship between spect ral types and 



effecti ve temperatures caused by non-solar metallicities, as for instance in the work by lKudritzki. Bresolin & Przyt 



(2003), are less important as long as the flux-weighted gravities log gf are well determined and the 
objects are hot enough. Errors are only introduced through erroneous bolometric corrections and 
errors in E(B-V), when the observed wavelength-dependent magnitudes have to be transformed 
into de-reddened bolometric magnitudes. While this may still be a serious error source and should 
be avoided through the application of our diagnostic technique described in the previous chapters, 
the uncertainties of the primary quantity in the application of the method are significantly reduced 
increasing the reliability of the FGLR method. 

Because of the importance of this effect for distance determinations we have investigated 
the physical reasons leading to Fig. [29] and [301 We discuss them in the remaining part of this 
subsection. The strength of absorption lines depends on the line strength /3 V 

j3 v = j4/<, (18) 

where k^ v is the line absorption coefficient and K v the continuum absorption coefficient at the 
frequency of the line. For H<5 and all other Balmer lines we have 



K"i v ocn 2 4>(Av) (19) 

with fi2 the occupation number in the second level of hydrogen and <!>( Av) the line broadening 
profile. Note that we have neglected the influence of stimulated emission, which is appropriate at 
the temperatures in question and for this qualitative discussion. For the supergiants the higher 
Balmer lines are always saturated in the Doppler line cores of the line profiles, where they simply 
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reflect the source function in the outer atmosphere. However, the contribution to the equivalent 
widths of these Doppler cores is small. More important is the contribution from the line wings, 
where the line broad ening is domin ated by th e Stark effect and the profile function is very well 
approximated by (see lUnsoeldl l 1968. p. 323 or [Mihalasll 19781 . p. 296) 



<t>(Av) n E AX~ 2 - 5 . (20) 

n E is the local electron density. The continuum absorption coefficient at the wavelength of 
is well described by 



K L v = n E o E (\+n E f(X,T)o/o E ). (21) 

a is a cross section of 10 -18 cm 2 and o E is the cross section for Thompson scattering at 
free electrons. The function f(X,T) contains the contributions of all continuous true absorption 
processes, but is dominated by hydrogen free-free and bound-free absorption. In LTE a good 
approximation (better than 10%) in the range between 8000 K and 15000 K and at the wavelength 
ofH5 is 



f{X J) « 1.06 Kr 19 r 4 -3 (22) 

in units of cm 3 . T4 is the temperature in units of 10 4 K. NLTE effects in the third level of 
hydrogen can have a small influence on f(X,T), but will change mostly the value of the constant 
but not the temperature dependence. The occupation number of hydrogen in the second level is 
given by 

n 2 ^b2n 2 E T- L5 exp(E 2 /kT) (23) 

with b 2 the non-LTE departure coefficient and E 2 the ionization energy from the second level, 
k is the Boltzmann constant. For this equation we have assumed that most of the free electrons in 
the atmospheres come from hydrogen, a valid approximation. In our temperature range we find 
that we can approximate (within 25%) 

T 4 L5 exp(E 2 /kT) » 6.310 l T 4 6 (24) 
This leads to the following expression for the line strength /3 V 



/3 V - (n E T 4 3 ) 2 /(\ + 1.06 10- 1 9 n E T 4 3 o/o E ). 



(25) 



Thus, the line strength in the wings of the hydrogen Balmer lines is obviously a function of 
electron density devided by the third power of temperature in the layers of the formation of line 
wings, which is close to the continuum optical depth x c v = 2/3. It is straightfoward to estimate 
n E T~ 3 in these layers. Starting from the hydrostatic equation 

dP/dx c v = -g(\-T)p/K L v (26) 

and with p = m p n p {\ +47), n p = n E , P = n E (2 + Y)kT and p, P, n p , Y being mass den- 
sity, gas pressure, proton number density and number ratio of helium to hydrogen, respectively, 
we obtain a differential equation for the electron density as a function of continuum optical depth 
(making the usual assumption that the logarithmic density gradient is much larger than the loga- 
rithmic temperature gradient) , which reads 

( 1 + n E f(X J)a/a E )dn E =A{g,TJ)dx c v: (27) 

where 

A(g, r, T) = g(l - T)m p (\ + 47)/(ofe(2 + Y)kT) (28) 

r is the usual Eddington-ratio of Thompson scattering radiative acceleration devided by grav- 
ity, namely (assuming again n E = n p ) 

r = 7.49r 4 4 (i+y)/(s(i+4y) -g- 1 ( 29) 

The integration of this differential equation is simple (again regarding the temperature gradi- 
ent as negligible). We obtain at every continuum optical depthr c 

n E + n 2 E af(X,T)/(2a E ) =A(gJJ)x c v (30) 

or 

n E of(Xj)/o E = (l + 2A(g,TJ)T c v of(Xj)/o E ) l l 2 -\. (31) 
With /(A, T) oc T 3 we see immediately that 



n E T^ 3 oc (\ +g F (\-T)* const. 



(32) 



which means that the line strength /3 V depends solely on the flux weighted gravity and not on 
effective temperature. We also note that depending on whether the second term under the square 
root is small or large compared with unity (i.e. gp(l — T) small or large) we will have heT 
£f(1 — r) or (#f(1 — r)) ' 5 , respectively. Numerical solutions of the exact problem including 
full NLTE opacities and the temperature gradient show that this simple estimate is qualitatively 
correct and fairly accurate. 

These simple calculations explain the physical reason why the the strengths of the Balmer 
lines depend mostly on the flux weighted gravity and why the isocontours of the equivalent widths 
in Fig. [30] are horizontal. However, they do not explain, why the isocontours bend upwards at 
higher temperature. The reason is the recombination of hydrogen at lower temperatures. While it 
is intuitively clear that recombination must lead to stronger hydrogen lines and thus to an upward 
turn of the isocontours, we verified this by simple numerical calculations solving the same problem 
as above but including the effects of recombination through the Saha formula. We also included the 
effects of atmospheric temperature stratifications by allowing for appropriate T(Ty) relationships. 
At lower T e g and higher log g hydrogen partly recombines. Ky becomes smaller (°c n\) and T v is 
reached in much larger geometrical depth of much higher pressure and total number density. As 
a result, the effect of recombination is overcompensated and ng at the optical depth of line wing 
formation is much larger than in the calculation without recombination at these temperatures and 
gravities. 



6.2. The FGLR of NGC 300 



With 24 A supergiants studied in this work and the 6 early type B supergiants investigated by 
Urbaneja et al.1 (I2005|) we have a large sample, which we can use for an observational test of the 
FGLR. For the A supergiants, all the data and their u ncertainties hav e been discussed in the previ- 
ous chapters. A simi lar discussion has be en given by lUrbaneja et al.l (|2005l) for the B supergiants. 



Note, however, that lUrbaneja et al.l (12005b used the ground-based photometry by iBresolin et al. 



(|2002l) . whic h as discussed above has been superseded by the HST/ACS photometry presented and 
discussed by IBresolin et al.1 (|2005f) . With the effecti ve temperatures , grav ities and intrinsic energy 
distributions determined for each B supergiant by lUrbaneja et al.l (|2005|) we have used the new 
photometry to re-determine interstellar reddening and ex tinction, which then y ields new stellar 
radii and absolute V and bolometric magnit udes. Note that lUrbaneja et al.1 (|2005r) also used the old 
distance modulus bv lFreedman et al. (2001), which is significantly different from the new distance 
modulus determined by iGieren et al.1 (|2005ai) based on on optical and IR observations of a large 
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sample of Cepheids minimizing the errors caused by reddening. 

The result is shown in Fig. [3TJ which reveals a clear and rather tight relationship of flux 
weighted gravity log gp with bolometric magnitude M^y. A simple linear regression yields bpcLR 
= 8.1 1 for the zero point and Q-fglr = -3.52 for the slope. The standard deviation from this rela- 
tionship is o = 0.34 mag. 

We also overplot the FGLRs predicted by stellar evolution theory ( Meynet & Maedeill2005 ) as 
discussed in the previous paragraphs. At lower luminosities (or large log gf) they agree well with 
the observed linear regression. In particular, the FGLR for Milky Way metallicity is very close. 
For higher luminosities (or small log gp) the linear regression and the stellar evolution FGLRs 
diverge because of the curvature of the latter, which is caused by the fact that the exponent of the 
mass-luminosity relationship becomes smaller. Looking at the scatter of the data points around 
either the linear regression or the FGLR for Milky Way metallicity it is hard to tell which of the 
two is the better fit. Approximating the stellar evolution FGLRs by 



M bol = a(log g F ) (log g F -xo) +M% ol 
a(log gp) = ocq + cci (xo - log g F ) + a 2 (x - log g F ) 2 



(33) 



(values for Oq, 0C\, a 2 and xq are given in Table 4) we can calculate standard deviations to 
characterize the scatter around those. For Milky Way metallicity we obtain o = 0.35 mag, practi- 
cally identical to scatter around the linear regression. The FGLR for SMC metallicity has a much 
stronger curvature (as the result of less mass-loss in the O-star stage of stellar evolution) and is 
shifted towards higher luminosities and seems to be worse fit of the data. This is confirmed by the 
standard deviation of o = 0.41 mag for this relationship. Applying a vertical shift of AM^ 0/ = 0.16 
mag to the SMC metallicity relationship minimizes its standard deviation to o = 0.38 mag, still 
worse than for the linear regression and the FGLR for Milky Way metallicity. We note that most 
of our objects have metallicities between the Milky Way and the SMC comparable to the LMC. 
Unfortunatly, no complete set of models including the effects of rotation in the mass range needed 
is available at this metallicity at the moment. 

The new linear regression coefficient s for th e FGLR in NGC 300 are different from those ob- 



tained by iKudritzki. Bresolin & Przybillal (I2003|) (cifglr = -3.85 and bpcLR = 7.96) resulting now 



in a higher luminosity zero point and a somewhat shallower slope. This is caused by a combination 
of different factors, a new distance modulus to NGC 300, improved HST photometry leading to 
larger extinction and a systematic change of stellar parameters as the result of the new diagnostic 
technique of the low resolution spectra. 
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6.3. The FGLR averaged over eight galaxies 



In their first investigation of the empirical FGLR lKudritzki. Bresolin & Przybillal (120031) have 
added A supergiants from six Local Group galaxies with stellar parameters obtained from quanti- 
tative studies of high resolution spectra (Milky Way, LMC, SMC, M31, M33, NGC 6822) to their 
results for NGC 300 to obtain a larger sample. They also added 4 objects from the spiral galaxy 
NGC 3621 (at 6.7 Mpc) which were studied at low resolution. We will add exactly the same data 
set to our new enlarged NGC 300 sam ple, however, with a fe w mi nor modifications. For the M ilky 
Way we include the latest results from lPrzybilla et al.l (|2006f) and lSchiller & Przv billa (2007) and 
for the two objects in M31 we use the new stellar parameters obtained by Przvbilla et al.1 (|2008 ). 
For the objects in NGC 3621 we apply new HST photometry. (Note that for the objects T e g and 
log g were obtained from a spectral fit assuming LMC metallicity, i.e. the detailed technique devel- 
oped here for the NGC 300 objects has not yet been applied. However, eliminating those objects 
from the sample does not change the results significantly). We have also re-analyzed the LMC 
objects using ionization equilibria for the temperature determination. The results are summarized 
in Table 5. The bolometric corrections for all objects are calculated using eq.(6). 



Fig. [32] shows bolometric magnitudes and flux- weighted gravities for this full sample of eight 
galaxies revealing a tight relationship over one order of magnitude in flux-weighted gravity. The 
linear regression coefficients are cifglr = —3.41 ±0.16 and bpGLR = 8.02 ±0.04, very simimilar 
to the NGC 300 sample alone. The standard deviation is o = 0.32 mag. The new coefficients are 
slightly different from the previous work by lKudritzki. Bresolin & Przybillal (|2003l) (cifglr = -3.71 
and bfGLR = 7.92) resulting in a FGLR brighter by 0.1 to 0.2 mag now in much closer agreement 
with stellar evolution. Compared with the new results the stellar evolution FGLR for Milky Way 
metallicity provides a fit of almost similar quality with a standard deviation of o = 0.31 mag. The 
stellar evolution FGLR fit for SMC metallicity is worse with o = 0.42 mag (applying AM^ 0/ = 0.20 
mag reduces the standard deviation to o = 0.35 mag). 

We conclude that the simple linear regression fit of this sample is the best way to describe 
the empirical FGLR at this point. It is basically in agreement with stellar evolution theory, at least 
if mass -loss and rotational mixing affect the evolution as for the models with solar metallicity. 
We note that for the evolutionary models accounting for mass-loss and mixing based on SMC 
metallicity we have an offset relative to the observed FGLR by about 0.20 mag and a curvature 
stronger than observed. Again, we stress the need for a complete set of models including the 
effects of rotation at LMC metallicities. 
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7. Conclusions and Future Work 

The goal of this work has been to demonstrate the astrophysical potential of low resolution 
spectroscopy of A supergiant stars in galaxies beyond the Local Group. By introducing a novel 
method for the quantitative spectral analysis we are able to determine accurate stellar parameters, 
which allow us to test stellar evolution models including the evolutionary time scales in the A su- 
pergiant stage. Through the spectroscopic determination of stellar parameters we can also constrain 
interstellar reddening and extinction by comparing the calculated SED with broad band photome- 
try. We find a very patchy extintion pattern as to be expected for a star forming spiral galaxy. The 
average extinction is in agreement with multi-wavelength studies of Cepheids including K-band 
photometry. 

The method also allows to determine stellar metallicities and to study stellar metallicity gradi- 
ents. We find a metallicity close to solar in the center of NGC 300 and a gradient of —0.08 dex/kpc. 
To our knowledge this is the first systematic stellar metallicity study in galaxies beyond the Local 
Group focussing on iron group elements. In the future the method can be extended to not only 
determine metallicity but also the ratio of a- to iron group elements as a function of galactocentric 
distance. The stellar metallicities obtained can be compared with oxygen abundance studies of 
HII regions using the strong line method. This allows us to discuss the various calibrations of the 
strong line method, which usually yield very different results. 

The improved spectral diagnostic method presented here enables us to very accurately de- 
termine stellar flux weighted gravities log gF = log g/Tjjf and bolometric magnitudes, which we 
explain by a detailed discussion of the physical background of the spectroscopic diagnostics. We 
find that above a certain threshold in effective temperature a simple measurement of the strengths 
of the Balmer lines can be used to determine accurate values of log gp. 

Absolute bolometric magnitudes Mt, i and flux-weighted gravities log gF are tightly corre- 
lated. It is shown that such a correlation is expected for stars, which evolve at constant luminosity 
and mass. We discuss the observed "flux- weighted gravity - luminosity relationship (FGLR)" in 
detail and compare with stellar evolution theory. We find resonable agreement with evolutionary 
tracks which account for mass-loss and rotational mixing assuming solar metallicity. The agree- 
ment is less good for similar tracks with SMC metallicity. 

With a relatively small residual scatter of o = 0.3 mag the observed FGLR is an excellent tool 
to determine accurate spectroscopic distance to galaxies. It requires multicolor photometry and 
low resolution (5A) spectroscopy to determine effective temperature and gravity and, thus, flux- 
weighed gravity directly from the spectrum using the self-consistent method developed here. With 
effective temperature, gravity and metallicity determined we also know the bolometric correction, 
which is small for A supergiants, which means that errors in the stellar parameters do not largely 
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affect the determination of bolometric magnitudes. Moreover, we know the intrinsic stellar SED 
and, therefore, can determine interstellar reddening and extinction from the multicolor photometry, 
which will allow the accurate determination of the reddening-free apparent bolometric magnitude. 
The application of the FGLR will then yield absolute magnitudes and, thus, the distance modulus. 
With the intrinsic scatter of a = 0.3 mag and 30 targets per galaxy one can estimate an accuracy of 
0.05 mag in distance modulus (0.1 mag for 10 target stars). 

The advantage of the FGLR method for distance determinations is its spectroscopic nature, 
which provides significantly more information about the physical status of the objects used for 
the distance determination than simple photometry methods. Most importantly, metallicity and 
interstellar extinction can be determined directly. The latter is crucial for spiral and irregular 
galaxies because of the intrinsic patchiness of reddening and extinction. 

Since supergiant stars are known to show intrinsic photometric variability, the question arises 
whether the FGLR method is affected by such variability. For the our targets in NGC 300 this 



issue has been carefully investigated by iBresolin et al.l (12004). who studied their CCD photometry 



lightcurves obtained over many epochs in the parallel search for Cepheids in NGC 300. They 
concluded that amplitudes of photometric variability are very small and do not affect distance 
determinations using the FGLR method. The standard deviations from the mean V magnitude 
for subset of A supergiants identified as variable ranges between 0.03 to 0.05 magnitudes, with 
maximum amplitudes of the variability between 0.08 to 0.23 mag. This is clearly within the one o 
uncertainty found for the FGLR. 

The effects of crowding and stellar multiplicity are also important. However, in this regard 
A supergiants offer tremendous a dvantages rela t ive to other stellar distance indicators. First of 



all, they are significantly brighter. IBresolin et al.l (|2005l) using HST ACS photometry compared to 



ground-based photometry have studied the effects of crowding on the Cepheid distance to NGC 
300 and concluded that they are negligible. With A supergiants being 3 to 6 magnitudes brighter 
than Cepheids it is clear that even with ground-based photometry only crowding is generally not 
an issue for these objects at the distance of NGC 300 and, of course, with HST photometry (and 
in the future JWST) one can reach much larger distances before crowding becomes important. 
In addition, any significant contribution by additional objects to the light of an A supergiant will 
become apparent in the spectrum, if the contaminators are not of a very similar spectral type, 
which is very unlikely because of the short evolutionary lifetime in the A supergiant stage. It 
is also important to note that A supergiants have evolutionary ages larger than 10 million years, 
which means that they have time to migrate into the field or that they are found in older clusters, 
which are usually less concentrated than the very young OB associations. 

It is evident that the type of work described in this paper can be in a straightforward way 
extended to the many spiral galaxies in the local volume at distances in the 4 to 7 Mpc range. 
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Pushing the method we estimate that with present day 8m to 10m class telescopes and the existing 
very efficient multi-object spectrographs one can reach down with sufficient S/N to V = 22.5 mag 
in two nights of observing time under very good conditions. For objects brighter than My = -8 
mag this means metallicities and distances can be determined out to distances of 12 Mpc (m-M = 
30.5 mag). This opens up a substantial volume of the local universe for metallicity and galactic 
evolution studies and independent distance determinations complementary to the existing methods. 
With the next generation of extremely large telescopes such as the TMT, GMT or the E-ELT the 
limiting magnitude can be pushed to V = 24.5 equivalent to distances of 30 Mpc (m-M = 32.5 
mag). 
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Fig. 1. — Parameter space of the model atmosphere grid used in this study. Different symbols 
refer to the microturbulence values used for the models. Open circles: v t = 8 km/s; open squares 
7 km/sec; open triangles 6 km/sec; filled squares 5 km/sec; filled circles 4 km/sec. At each grid 
point metallicities from between [Z] = 0.3 to -1.30 are available (see text). 
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Fig. 2. — The temperature-metallicity degeneracy of spectral types. Synthetic spectra of three A 
supergiant model atmospheres degraded to 5 A resolution with {T e g, log g, [Z]) = (10500 K, 1.40, 
0.15) (top); (9500 K, 1.20, -0.30) (middle); (8750 K, 1.00, -0.50) (bottom). The latter two are 
vertically shifted by 0.1 and 0.2. Metal line identifications are given by the bars at the top (solid: 
Fell; dashed: CrII; dotted: Till; dashed-dotted: Mgll). Note the change of the strength of the Hel 
line at 4471 A from the hottest to the cooler models. 




Fig. 3. — Energy distributions of the models of Fig. [2} The fluxes are normalized at 4380 A. The 
hottest model has the steepest SED. 
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Fig. 4. — Balmer jump fluxes of the models of Fig. [2] The fluxes are normalized at 3790 A. The 
hottest model has the smallest Balmer jump. 
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Fig. 5. — Isocontours of H8 (solid) and Dg (dashed) in the {log gjog T e g) plane. H8 isocontours 
start with 1 A equivalent width and increase in steps of 0.5 A. Dg isocontours start with 0. 1 dex 
and increase by 0. 1 dex. The iscontours are calculated fro the model atmosphere and non-LTE line 
formation grid described in section 3 and Fig. 1. 
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Fig. 6. — Model atmosphere fit of two observed Balmer lines of target No. 21 of Table 1 for T e g 
= 10000 K and log g = 1.55 (solid). Two additional models with same T eff but log g = 1.45 and 
1.65, respectively, are also shown (dashed). Fits of other Balmer lines such as H r .8,9,io are similar. 
Note that the calculated line profiles were folded with a Gaussian of 5 A FWHM to account for 
instrumental broadening 
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Fig. 7. — Model atmosphere fit of the observed Balmer jump of target No. 21 of Table 1 for T e g 
= 10000 K and log g = 1.55 (solid). Two additional models with the same log g but T e g = 9750 K 
(dashed) and 10500 K (dotted) are also shown. The horizontal bar at 3600 A represents the average 
of the flux logarithm over this wavelength interval, which is used to measure D B (see text). 




Fig. 8. — (log g.log T e ff)-ht diagram for target No. 21 of Table 1. The solid curves represent the 
best fits of the Balmer jump and the Balmer lines. The dashed curves correspond to the maximum 
fitting errors. 




Fig. 9. — Same as Fig.[8]but for the flux weighted gravity log gF instead of gravity log g. 
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Fig. 10. — Effective temperature uncertainties ±AT e ff/T e ff as a function of effective temperatures 
(see Table 1). Open and solid symbols refer to positive and negative errors, respectively. 



-44- 




Fig. 11. — Synthetic metal line spectra calculated for the stellar parameters of target No. 21 as a 
function of metallicity in the spectral window from 4497 A to 4607 A. Metallicities range from 
[Z] = -1.30 to 0.30, as described in the text. The dashed vertical lines give the edges of the 
spectral window as used for a determination of metallicity. The synthetic spectra were folded 
with a Gaussian of 5 A FWHM to account for the instrumental broadening. Rotational broadening 
is negligible for A supergiants at this resolution and not taken into account. 
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Fig. 12. — Observed spectrum of target No. 21 for the same spectral window as Fig-fTTIoverplot- 
ted by the same synthetic spectra for each metallicity separately. The line identification symbols 
plotted for [Z] =0.15 are the same as in Fig. [TO For the renormalization of the observed spectrum 
at each different metallicity see discussion in the text. The abscissae are wavelengths in A. 
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Fig. 13. — X([Z])i as obtained from the comparison of observed and calculated spectra in Fig. IT2l 
The solid curve is a third order polynomial fit. 
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Fig. 14. — Determination of metallicity for target No. 21 in six more spectral windows. Line 
identifications are again given for the plots with [Z] =0.15. Solid bars refer to Fell, dotted bars 
represent Till in windows 2, 3, 4, 6, Sill in window 1 and Hel in window 5. Dashed bars indicate 
Hel in window 1 and 3 and Mgll in window 4 and 6. The dashed-dotted bar in window 4 is Hel. 
The abscissae are wavelengths in A. 
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Fig. 15. — X([Z])i for all seven spectral windows of target No.21. 
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Fig. 17. — Effective temperatures according to spectral types compared with effective temperatures 
obtained from the quantitative spectral analysis. For discussion see text. 
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Fig. 18. — Relative difference of spectral type and spectral analysis temperatures as a function of 
effective temperatures obtained from the quantitative spectral analysis. 
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Fig. 19. — Relative difference of spectral type and analysis temperatures as a function of metallicity 
obtained from the quantitative spectral analysis. 
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Fig. 20. — NGC 30 A supergiants (circles, this study) and early B supergiants (squares, 



Urbaneia et al.1 120051) in the (log g,log T e g) plane compared with evolutionary tracks by 
Mevnet & Maederl d20Q5h of stars with 15 M Q (solid), 25 M (dashed), and 40 M Q (dashed- 



dotted), r espectively. The tracks in clude the effects of rotation and are calculated for SMC metal- 
licity (see uVIeynet & Maederj|2005|) . 
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Fig. 21. — NGC 300 A supergiants (circles, this study) and early B supergiants (squares, 



Urbaneja et al.ll2005h in the HRD compared with evolutionary tracks for stars with 15 M (solid), 



20 M (long-dashed), 25 M Q (short-dashed), and 40 M Q (dashed-dotted) , respectively. The track s 
include the effects of rotation and are calculated for SMC metallicity (see Meynet & Maedej2005 ). 
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Fig. 22. — Comparison of evolutionary with spectroscopic masses for the NGC 300 A supergiants 
of this study 
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Fig. 23. — Logarithm of the ratio of evolutionary to spectroscopic masses as a function of the 
logarithm of spectroscopic mass. 
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Fig. 24. — Metallicity [Z] as a function of angular galacto-centric distance p/po ( see Table 1) for 



the A- supergiants of this work (filled circles) and the early B-supergiants studied by lUrbaneja et al. 



(|2005l) (filled squares). Note that for the latter metallicity refers to oxygen only. The dashed curve 
represents the regression discussed in the text. Note that typical uncertainties in [Z] are ±0.2 dex. 
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Fig. 25. — HST/ACS color magnitude diagram of the observed fi elds in NGC 300. The p osition of 
our targets is indicated by filled circles and evolutionary tracks (jMeynet & Maedern2005r) with 12, 
15, 20 and 25 M are overplotted. 
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Fig. 26. — The number of stars in Fig. [25] in 0.1 magnitude bins as a function of magnitude. The 
fit discussed in the text is shown as a straight line. 
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Fig. 27. — The FGLR as predicted by stellar evolution theory (IMeynet & Maederll2005|) using the 
mass-luminosity relationship of eq. 8 and the coefficients of Table 4. The solid curve is calculated 
for solar metallicity and includes the effects of stellar rotation. The dashed-dotted curve is for the 
same metallicity but does not include rotation. The long-dashed curves assumes SMC metallicity 
and accounts for rotation (for this metallicity the curve without rotation is very similar). 
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Fig. 28. — Equivalent width of H5 as a function of gravity log g for models of different effective 
temperature ranging between 8300 K (top curve, crosses) to 15000 K (bottom curve, filled circles). 
The Tgff steps can be inferred from Fig.Q] The metallicity is [Z] = -0.5, but the curves at different 
metallicity are very similar. 




Fig. 29.— 



Same as Fig. [28l but for the flux weighted gravity log gF- 




Fig. 30. — Same as Fig.[5]but for the flux weighted gravity log gF instead of gravity log g. 
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Fig. 31. — The FGLR of A (solid circles) and B (open circles) supergiants in NGC 300 and the 
linear regression (solid). The stellar evolution FGLRs of Fig. |27]for models with rotation are also 
overplotted (dashed: Milky Way metallicity, long-dashed: SMC metallicity). 
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Fig. 32. — Same as Fig.[3H but with additional objects from 7 other galaxies (see text). 
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Table 1 . NGC 300 - Spectroscopic targets and stellar parameters 



No. 


name 


P/Po 

a 


spec, 
type 


D B 


e J j 
K 


log g 

CgS 


log gf 
CgS 


[Z] 


E(B-V) 
mag 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 





A8 


1.03 


B9-A0 




loooo+S!! 

— 3UU 


1.60+°" 


1.60+ , 08 


-0.44±0.15 


0.09 


b 


1 


A10 


0.81 


A2 


0.37 


9250+?™ 


1 45+0.13 


1 59+o° 8 

—0.06 


-0.04±0.20 


0.16 




2 


All 


0.94 


B8 


0.08 


10500+^00 


1-251" 


2 17+0.05 


-0.65±0.15 


0.16 


d 


3 


C6 


0.85 


AO 


0.34 


9750+^0 

—400 


1.70+°" 

—0.18 


1.74+°+° 

—0.10 


-0.45±0.20 


0.16 




4 


C16 


0.70 


B9 


0.13 


10500+500 

— DUU 


1.35+°' 2 


1.27+0-04 

— U.vj 


-0.43±0.15 


0.14 


cj 


5 


D12 


0.64 


A2 


0.28 


9000+4°,°, 


1.10l° : J7 


1.28+°+° 


-0.16±0.15 


0.24 




6 


A6 


0.98 


A1-A2 


0.56 


8500+J™ 


1.50+° : 22 


1.78+°+ 7 


-0.48±0.30 


0.17 




7 


A13 


0.83 


B8 


0.22 


11500+™- 


1 95+° 17 


i 71+O.O6 
1 -' 1 -0.07 


-0.48±0.30 


0.11 




8 


A18 


0.93 


B8 


0.19 


12000^0 


2.oo + °;!^ 


1 fiR+0' 06 
L(> 8-0.06 


-0.36±0.30 


0.13 


C 


9 


B5 


0.35 


B8 


0.23 


11000+5™ 


1.80+°;J 4 


1 63+° 05 

1 - DJ -0.05 


-0.29±0.20 


0.11 


c 


10 


B8 


0.26 


AO 




9 500+ 400 


1 65+° 19 


i 74+0.12 
1 -'+-0.13 


-0.14±0.25 


0.12 


b,c 


11 


BIO 


0.14 


A2-A3 


0.40 


9250+ 530 


1 60+° 32 

l.OU_ 24 


, 7S +0.20 
1' 78 -0.15 


-0.07±0.30 


0.15 




12 


Bll 


0.13 


A4 


0.54 


8500+J3 , 


1 50+° 27 


1 74+0.22 


-0.09±0.30 


0.11 




13 


B19 


0.36 


B8-B9 


0.26 


10500+|° 


, 7C+O.I6 
1 -'- 3 -0.18 


1 67+° 08 
1 -°'-0.08 


-0.30±0.20 


0.15 




14 


CI 


1.18 


Al 


0.49 


8750t!?° 


1 50+°' 19 


, 7 -!+0.14 
l. ; J_ .14 


-0.33±0.25 


0.15 




15 


C8 


0.92 


B9 


0.19 


H500+™ 




1 61+ 006 

1 - oi -0.07 


-0.45±0.30 


0.12 




16 


C9 


0.93 


B9 


0.39 


9500+«° 


1 70+ a21 
1 - /u -0.20 


1 7Q+0.13 
1 - /y -0.13 


-0.44±0.30 


0.12 




17 


C12 


0.90 


B9-A0 


0.49 


8750+J™ 


1 45+0.18 


1 fiS+013 
1(> 8-0.16 


-0.54±0.25 


0.09 




18 


D2 


0.57 


Al 


0.41 


9250+400 


1 60+ a21 


i 74+0.13 


-0.41±0.20 


0.16 




19 


D7 


0.71 


A4 


0.56 


8500+1™ 


1 70+0' 24 
1 - /u -0.20 


1 98+° 19 
1 - yo -0.16 


-0.31±0.20 


0.23 


c 


20 


D10 


0.57 


A4 


0.63 


8300+270 


1 60+ a26 
1 - DU -0.26 


I 97+0 21 
1 - yz -0.21 


-0.27±0.15 


0.18 


c 


21 


D13 


0.79 


AO 


0.27 


10000+j°° 


1 55+ al9 


, cc+0.09 
1 - JJ -0.08 


-0.39±0.15 


0.13 




22 


D17 


0.82 


B9-A0 


0.32 


9750+J™ 


1 60+°' 12 


1 64+°' 08 
^ -0.08 


-0.46±0.15 


0.15 




23 


D18 


1.08 


B9 


0.41 


9250+*°, 


1 60+°+ 8 
lou — 0.13 


i 74+0.12 


-0.66±0.20 


0.07 





"Galactocentric distance, in units of the isophotal radius po = 9.75 arcmin ~ 5.33 kpc 
b no Balmer jump measured; T e tr from Hel lines 

c no HST photometry available; ground-based photometry used, see text 
d HeI lines also used for used T e ff, see text 
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Table 2. Bolometric corrections and magnitudes, radii and masses 



No. 


name 


in y 


My 


BC 


M OQ \ 


R/R 


M/M:, 


M/M 






mag 


" mag 


mag 


a mag 




spec 


evol 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 





A8 


19.41 


-7.24 


-0.29 


-7.53±0.13 


96 


13.2 


15.7 


1 


A10 


18.93 


-7.94 


-0.12 


-8.06±0.05 


142 


20.8 


18.8 


2 


All 


18.27 


-8.59 


-0.49 


-9.08±0.12 


177 


20.2 


27.7 


3 


C6 


19.84 


-7.01 


-0.26 


-7.27±0.12 


89 


14.4 


14.4 


4 


C16 


18.07 


-8.75 


-0.44 


-9.19±0.13 


186 


28.1 


29.0 


5 


D12 


18.66 


-8.45 


-0.13 


-8.58±0.04 


191 


16.7 


22.7 


6 


A6 


19.78 


-7.11 


-0.03 


-7.14±0.05 


110 


14.0 


13.9 


7 


A13 


19.85 


-6.86 


-0.61 


-7.47±0.13 


70 


16.0 


15.4 


8 


A18 


19.99 


-6.79 


-0.66 


-7.45±0.17 


64 


14.8 


15.3 


9 


B5 


19.77 


-6.94 


-0.48 


-7.42±0.11 


75 


12.9 


15.2 


10 


B8 


19.83 


-6.89 


-0.19 


-7.08±0.12 


86 


12.0 


13.6 


11 


BIO 


19.69 


-7.14 


-0.14 


-7.28±0.10 


99 


14.3 


14.5 


12 


BI1 


19.95 


-6.76 


0.00 


-6.76±0.06 


93 


9.9 


12.3 


13 


B19 


20.23 


-6.60 


-0.38 


-6.98±0.06 


67 


9.2 


13.2 


14 


CI 


18.99 


-7.85 


-0.06 


-7.91±0.05 


148 


25.3 


17.8 


15 


C8 


20.01 


-6.74 


-0.57 


-7.31±0.16 


65 


10.9 


14.6 


16 


C9 


20.22 


-6.52 


-0.21 


-6.73±0.09 


73 


9.7 


12.2 


17 


C12 


20.34 


-6.31 


-0.06 


-6.37±0.04 


73 


5.5 


11.0 


18 


D2 


19.83 


-7.04 


-0.16 


-7.20±0.07 


96 


13.3 


14.1 


19 


D7 


20.02 


-7.06 


-0.01 


-7.07±0.05 


107 


20.8 


13.6 


20 


D10 


19.75 


-7.18 


0.02 


-7.16zb0.05 


117 


19.7 


13.9 


21 


D13 


18.86 


-7.92 


-0.29 


-8.21±0.11 


131 


22.0 


19.8 


22 


D17 


19.44 


-7.40 


-0.25 


-7.65zb0.05 


106 


16.3 


16.3 


23 


D18 


19.96 


-6.62 


-0.18 


-6.80±0.06 


80 


9.2 


12.5 



"distance modulus = 26.37 mag, iGieren et al] l2005d) 
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Table 3. Central metallicity [Z] and metallicity gradient in NGC 300 



Source 


Central abundance 


Gradient 


Gradient 


Comments 




dex 


P/Po 


dex/kpc 




(1) 


(2) 


(3) 


(4) 


(5) 


Dopita & Evans(1986) 


0.26±0.17 


-0.63±0.10 


-0.118±0.019 


HII, oxygen 


Zaritsky et al. (1994) 


0.29±0.04 


-0.54±0.09 


-0.101±0.017 


HII, oxygen 


Kobulnicky et al. (1999) 


0.06±0.04 


-0.27±0.09 


-0.051±0.017 


HII, oxygen 


Denicolo et al. (2002) 


-0.08±0.05 


-0.46±0.10 


-0.086±0.019 


HII, oxygen 


Pilyugin (2001) 


-0.17±0.06 


-0.28±0.12 


-0.053±0.023 


HII, oxygen 


Pettini & Pagel (2004) 


-0.19±0.04 


-0.36±0.08 


-0.068±0.015 


HII, oxygen 


This work 


-0.06±0.09 


-0.44±0.12 


-0.083±0.023 


stars, metals 



Table 4. Fit coefficients for stellar evolution mass -luminosity relationship and FGLR 



Model 


Metallicity 


a 


b 


c 


log(M min /M (S ) 


ao 




a 2 


xo 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


with rotation 


Milky Way 


3.80 


1.40 


4.095 


0.943 


3.505 


0.15 


0.30 


1.957 


-6.464 


with rotation 


SMC 


3.10 


1.20 


4.560 


1.077 


3.659 


0.32 


0.32 


1.967 


-6.439 


no rotation 


Milky Way 


4.25 


2.00 


3.825 


0.943 












no rotation 


SMC 


4.10 


2.05 


4.030 


0.954 
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Table 5. Stellar parameters of A supergiants in 7 galaxies 



No. 


name 


galaxy 


m-M 


T eff 


logg 


Mbol 


comments 








mag 


K 


cgs 


mag 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 





Tj Leo 


Milky Way 


9.00 


9600 


2.00 


-5.73 




1 


riDl 1 lol j 


MilKy Way 


1 1.80 


o i cn 
9150 


1.45 


-7.43 




1 


riDyzzU / 


MilKy Way 


12.40 


9500 


1.20 


O A/1 

-9.04 




•2 
J 


Rigel 


MilKy Way 


7.80 


1 2000 


1.75 


-8.49 


a 


1 

4 


Deneb 


Milky Way 


O 

9.50 


8525 


1.10 


-8.48 


b 


c 

J 


HUlzyjo 


MilKy Way 


1 1.80 


9250 


1.20 


-8.08 


c,d 


J 


rlL>144oy 


MilKy Way 


1 1.80 


9000 


1 .35 


-7.81 


c,d 


5 


HDZZjjoj 


MilKy Way 


12.00 


8500 


0.95 


-5.43 


c,d 





oi. /;n oil 
oK-oy z 1 1 


LMC 


18.50 


1 1000 


1.30 


n an 
-9.00 


e 


7 


III 


LML 


18.50 


10500 


1.20 


-9.53 


e 


8 


Sk-67 207 


LMC 


18.50 


10000 


1.20 


-8.50 


e 


9 


Sk-66 58 


LMC 


18.50 


9500 


1.10 


-8.74 


e 


10 


Sk-67 201 


LMC 


18.50 


10000 


1.20 


-9.14 


6 


11 


Sk-69 239 


LMC 


18.50 


10000 


1.05 


-9.68 


e 


12 


Sk-69 170 


LMC 


18.50 


10000 


1.20 


-8.89 


e 


13 


Sk-69 299 


LMC 


18.50 


8750 


1.00 


-8.93 


e 


14 


Sk-67 44 


LMC 


18.50 


8000 


0.65 


-9.54 


e 


15 


AV457 


SMC 


19.00 


8500 


0.90 


-9.10 


e 


16 


41-3654 


M31 


24.47 


9200 


1.00 


-9.08 


r.k 


16 


41-3712 


M31 


24.47 


8850 


1.00 


-8.79 


r.k 


17 


11 7- A 


M33 


24.47 


9500 


1.10 


-8.29 


%4 


18 


NGC6822m 


NGC 6822 


23.31 


9000 


1.40 


-7.23 


h.i 


19 


No. 1 


NGC 3621 


29.08 


9250 


1.30 


-8.21 


&,j 


20 


No. 9 


NGC 3621 


29.08 


9250 


1.12 


-9.20 


dj 


21 


No. 16 


NGC 3621 


29.08 


10500 


1.55 


-8.05 


i,j 


22 


No. 17 


NGC 3621 


29.08 


8350 


1.15 


-8.58 





tezvbilla etai]j2006l) 

ISchiller & Przvbilla I UtKfh 

■ kudritzki et~aD jl999h 

iKudritzki. Bresolin & Przvbilhj JiooH) 

e this work 

1Przvbillaetai]j2008h 

g distance modulus from lBonanos et aTl 120061) 

h distance modulus from Gier en et all 120061) 

IPrzvbillaliiobl) 

jBresolinetai]|200ll) 

k distance modulus from Holland] Il998t) 



M_bol 




M_bol 




